NASA— TM— 111597 





AIAA 96-1800 

Analytic Corrections to CFD Heating 
Predictions Accounting for Changes 
in Surface Catalysis 

Peter A. Gnoffo 

NASA Langley Research Center 
Hampton, VA 23681-0001 

George R. Inger 
Iowa State University 
Ames, LA 50011-3231 


31st AIAA Thermophysics Conference 
June 17-20, 1996/New Orleans, LA 


For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 
370 L'Enfant Promenade, S.W., Washington, D.C. 20024 




Analytic Corrections to CFD Heating Predictions Accounting for Changes in Surface Catalysis 


Peter A. Gnoffo * 

NASA Langley Research Center 
Hampton, Virginia 23681-0001 
and 

George R. Inger t 
Iowa State University 
Ames, Iowa 50011-3231 


Abstract 


Integral boundary-layer solution techniques ap- 
plicable to the problem of determining aerodynamic 
heating rates of hypersonic vehicles in the vicinity of 
stagnation points and windward centerlines are briefly 
summarized. A new approach for combining the in- 
sight afforded by integral boundary-layer analysis with 
comprehensive (but time intensive) computational fluid 
dynamic (CFD) flowfield solutions of the thin-layer 
Navier-Stokes equations is described. The approach 
extracts CFD derived quantities at the wall and at the 
boundary layer edge for inclusion in a post-processing 
boundary-layer analysis. It allows a designer at a work- 
station to address two questions, given a single CFD 
solution. (1) How much does the heating change for a 
thermal protection system with different catalytic prop- 
erties than was used in the original CFD solution? (2) 
How does the heating change at the interface of two dif- 
ferent TPS materials with an abrupt change in catalytic 
efficiency? The answer to the second question is partic- 
ularly important, because abrupt changes from low to 
high catalytic efficiency can lead to localized increase 
in heating which exceeds the usually conservative esti- 
mate provided by a fully catalytic wall assumption. 
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Nomenclature 

C relaxation coefficient in Method 2 

Cd equilibrium constant, p e e° D 

C p specific heat of molecules, J/kmoie-K 
h enthalpy per unit mass of mixture, J/kg 

ho heat of formation per unit mass of species, J/kg 

i(i-3,D) reaction integrals 

k r recombination rate coefficient, m 6 /kmole 2 -s 

Le Lewis number 

p pressure, N/m 2 

Pr Prandtl number 

q surface heating, W/m 2 

Q heat transfer function 

r ratio of velocity gradients 

R radius of curvature, m 

R s gas constant for species s J/kmole-K 
Ru univ. gas constant, 8314.3 J/kmole-K 

Sc Schmidt number 

T temperature, K 

Td characteristic temperature of dissociation, K 

it,w streamwise and crossflow velocites, m/s 

V total velocity, m/s 

x, z streamwise and crossflow distances, m 

Z bridging function 

a mass fraction of atoms 

P velocity gradient, 1 /s 

7 catalytic efficiency of surface 

r* composite Dahmkohler number 

T c , f c Dahmkohler numbers for surface catalysis 

Td Dahmkohler number for dissociation dominated flow 

Tq Dahmkohler number for recombination dominated flow 

p density, kg/m 3 

(i viscosity, N-s/m 2 

0 temperature ratio T/T e 

6d temperature ratio To fT e 

r} Levy-Lees transformation coordinates 

X parameter u e (d£/dx)( 1 -f r)/£, 1 /s 
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Superscripts 

+ value immediately following jump 

— value immediately preceding jump 


Subscripts 


oo 

free stream 

0 

computed using 7 = 0 

1 

computed using 7 = 1 

e 

boundary-layer edge 

EQ 

equilibrium 

F 

frozen 

hi 

CFD mesh point location 

J 

location of jump in surface catalysis 

N 

nitrogen atom 

O 

oxygen atom 

w 

wall 

7 

computed using 7 = 7 (T) 

7 1 

computed across discontinuity in 7 
from 7 = 7 (T) to 7 = 1 


Introduction 

Integral boundary-layer solution techniques ap- 
plicable to the problem of determining aerodynamic 
heating rates of hypersonic vehicles in the vicinity of 
stagnation points, windward centerlines, and swept- 
wing leading edges are discussed in the literature . 1-4 
The analyses include effects of finite-rate gas chemistry 
across the boundary layer and finite-rate catalysis of 
atom recombination at the surface. A new approach for 
combining the insight afforded by integral boundary- 
layer analysis with comprehensive (but time intensive) 
computational fluid dynamic (CFD) flowfield solutions 
of the thin-layer Navier-Stokes equations is described. 
The approach extracts CFD derived quantities at the 
wall and at the boundary layer edge for inclusion in 
a post-processing boundary-layer analysis. The post- 
processed data base allows a designer at a workstation 
to address two questions, given a single CFD solution. 
(1) How much does the heating change for a thermal 
protection system with different catalytic properties 
then was used in the original CFD solution? (2) How 
does the heating change at the interface of two differ- 
ent TPS materials with an abrupt change in catalytic 
efficiency? The answer to the second question is partic- 
ularly important, because abrupt changes from low to 
high catalytic efficiency can lead to localized increase 
in heating which exceeds the usually conservative es- 
timate provided by a fully catalytic wall assumption. 
For a given trajectory point, the approach uses a single 


CFD solution obtained with a known variation of cat- 
alytic efficiency over the vehicle surface. Changes to 
CFD baseline heating levels are calculated as a func- 
tion of changes in catalytic efficiency derived from in- 
tegral boundary-layer solution techniques that utilize 
CFD generated edge and wall conditions. 

The present paper focuses on flow over a sphere at 
several trajectory points characteristic of reentry of a 
representative winged Single Stage to Orbit (SSTO) ve- 
hicle. CFD calculations are made for three wall catal- 
ysis models at each of these points. The models in- 
clude noncatalytic, finite catalytic (Shuttle tile), and 
fully catalytic wall boundary conditions. The post pro- 
cessed integral boundary- layer corrections using each 
of three baseline CFD heating results are compared 
to actual CFD calculations for the corresponding off- 
baseline case. 

Also presented are the results of preliminary tests 
on vehicles of realistic geometric complexity appropri- 
ate for Access-to-Space studies. Relative differences in 
predicted heating levels between single CFD runs with 
post-processed corrections and multiple CFD runs us- 
ing different wall catalysis models are discussed. 

The goal of this research is to create a stand-alone 
post- processing tool which can be used in preliminary 
design of thermal protection systems (TPS) for hyper- 
sonic vehicles. The tool would make extensive use of 
a small number of CFD solutions computed using a 
relatively simple surface catalysis model. Design iter- 
ations are conducted without need of additional CFD 
runs until convergence on a single concept is achieved, 
at which point CFD could be used to provide a final 
check and/or recalibration point. 


Baseline CFD Solution Algorithm 

The computational fluid dynamic (CFD) baseline 
solutions for real gas, viscous analysis are provided by 
the Langley Aerothermodynamic Upwind Relaxation 
Algorithm (LAURA ). 5-7 Comparisons to experimen- 
tal data for hypersonic flows in air are documented in 
the literature . 8-12 The code employs upwind-biased, 
point-implicit relaxation. Inviscid fluxes are approxi- 
mated with Roe’s averaging , 13 eigenvalue limiting (sim- 
ilar to Harten 14 ), and Yee’s symmetric total variation 
diminishing scheme . 15 Viscous fluxes are approximated 
with central differences. A model for surface catalytic 
efficiency 16 j(T) used in the present work is defined by 

70 = 40.e- 1144O -/ T - 1435 < T w < 1580 

7 o = 39.9 iO- 9 e 2141 °V^ 1580 <T W < 1845 

( 1 ) 

7 N = 0.061e~ 2480 1410 < T w < 1640 

77V = 0.00061e 509 ° / T - 1640 <T W < 1905 
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The catalytic recombination rate I\ w s for species s is 
then defined by 

K Wt , = 7s(T w WR s T w /27r (2) 


Integral Boundary- Layer Method 


Highlights 

Integral bound ary- layer theory for evaluating 
nonequilibrium effects on surface heating has been de- 
scribed previously. 2 ^ Key features of the analysis are 
reviewed here. First, the model accounts for both fi- 
nite catalysis of the surface and finite reactivity of the 
boundary layer. The finite catalysis of the surface en- 
ters the analysis through the Dahmkohler number, F c , 
the ratio of atom recombination time at the surface and 
a characteristic local diffusion time. 


r c 


2 p u 


1 1/2 


PwX 


ScK v 


( 3 ) 


The finite reactivity of the boundary layer enters 
the analysis through the Dahmkohler number, Tg, and 
T# the ratio of characteristic local flow time to gas 
phase reaction time. In the case of recombination dom- 
inated boundary layers (near stagnation points) 


Tg 


4 krTrM 

Rlx 


( 4 ) 


In the case of dissociation dominated boundary lay- 
ers (near windside centerline where viscous shear raises 
temperatures in the boundary layer) 


Td = — Tg (5) 

Pe 

Second, analytic solutions to the governing 
boundary-layer equations can be made with simplify- 
ing approximations in the vicinity of the vehicle stag- 
nation point (nose region), swept leading edge region, 
and windward centerline region. The parameter x, 


x = « e (dC/dx)(l + r)/£ (6) 

is a function of the Levy-Lees transformed coordinate £ 
and takes on relatively simple, limiting values in these 
regions. The Levy-Lees transformed coordinate is de- 
fined by 

£ = 

Three-dimensional effects are included with the param- 
eter 1 + r, where r is the ratio of velocity gradients at 
the edge of the boundary layer. 

r = (dw e /dz)/(du e /dx) (8) 


/ 


pu) Pw ^ e dx . 


( 7 ) 


Global Change in Nose Region - Method 1 

The post-processing algorithm used to correct for 
changes in surface heating associated with global 
change of the surface catalysis model in the nose re- 
gion of a hypersonic vehicle is denoted Method 1. Sur- 
face heating in the nose region can be expressed 2 as an 
appropriate interpolation of an equilibrium boundary- 
layer heating rate Q Wj eq and a frozen boundary-layer 
heating rate Q Wy r by 

Qw — Qw,eq + Z(r)(Q WlF -Q WtEQ ) (9) 

In Eq. 9 the variable Q is a dimensionless heating rate, 
related to the dimensional heating rate q w by 


Qw 


— Pr/CpT e 

q w 

VPuJhvX 


(10) 


The limiting form of Eq. 10 for the stagnation point 
can be expressed by 


Qw ,i ,j 


Qw,i,j 


— Pr w ,i J /CpT e , i J 

Pw,i ,j Pw f i,jPe,iJ (1 P ij ) 


( 11 ) 


where Pe,i,j{ 1+r) is the limiting form of x in the stagna- 
tion region and 8 e ,i,j is the streamwise velocity gradient 
which, using Newtonian theory, can be approximated 
by 


Pi 


*>*,j 


du et j 

dx 


R. 


■r,tj 


f 2(Pe T i,jf Poo) 

Pe,i,j 


( 12 ) 


While Eqs. 11 and 12 are rigorously true only at the 
stagnation point, they are used throughout the nose 
region at every mesh point (i, j) of the baseline CFD 
solution in subsequent development of Method 1 . 

The equilibrium boundary-layer heating rate 
Qw,eq reference value at each mesh point (z, j) of the 
nose region is defined 2 by 

Qw,eq, ij = OA7PrH?j(l-0 w 

+ < 13 > 

The reference value for frozen bound ary- layer heating 
rate Q Wi f at each mesh point (i, j) of the nose region 
is defined 17 by 


Q w ,F,i,j = VAlPr'J* . (1 - e w 

*ehr> f 9 \ 
C p T e 1+r cJij 


+VPr x 


u c 4 . r e 0.67 a* 

’ 2 C p T e ^ C, 


(14) 


For stagnation region flow, the normalized kinetic en- 

u 2 ■ ■ 

ergy term * is small and can be omitted from 

both Eq. 13 and 14. Equation 14 accounts for fi- 
nite surface catalysis through the Dahmkohler number 
f ci j = T CyiJ /(0A7Sc 1 ^) but does not include effects 
of reactions across the boundary layer. 
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The bridging function Z(F*) in Eq. 9 varies from 0 
(equilibrium limit with three body recombination rates 
much faster than local flow rates) to 1 (frozen limit with 
three body recombination rates much slower than local 
flow rates). It is derived from a fundamentally based 
analysis of the atomic species conservation equations 
and is given by 


z(n 


3^1 + I6r* ./9 - 1 

2 T 4 rr 


( 15 ) 


where rj^ is a composite Dahmkohler number defined 
at each surface mesh point by 


r* = 


a 




h + 2 f eh + f 
(l + f c ) 2 


r G,ij (i6) 


■ *,j 


The integrals Ii, I 2 , and I 3 account for integrated ef- 
fects of reactions across the boundary layer on heating 
rates. They are tabulated 3 and curve fits are provided. 


h * 4 . 80 ( 0 . 50 / 5 ^)° 45 ^- 80(1 - w) - 1 84 

h « 1 . 80 ( 0 . 50 / 5 c u; ) o ' 12 ^°- 63(1 -“ ) - 1 15 ( 17 ) 

I 3 « O . 93 ( O . 5 O / 5 cu ,) - o ' 22 0 ° 41 ( 1-u '^“°- 65 


These integral approximations are valid over a param- 
eter range 0.04 < 9 W < 0.50 and — 2 < w < 0. In the 
test problems considered here it has been assumed that 
recombination of atomic oxygen can be used to char- 
acterize the reaction integrals in Eq. 17; consequently, 
k r = b.fllO 11 !^ with u> = —1.0. Values of the bridg- 
ing function in Eq. 15 were found to be substantially 
characteristic of frozen boundary- layer flow for all four 
points in Table 1. 

The following approximations have been employed 
in Method 1 to simplify extraction of CFD derived 
quantities. Because energy is used as a primary vari- 
able and temperature is a derived quantity in the CFD 
analysis we substitute h — ahn for C P T where ahjr, = 
and we substitute ( h — a/ip)/(/i — aho) e 
for 0 in the above equations. In order to model the 
catalytic effect of the surface on homogeneous recom- 
bimation of both nitrogen and oxygen atoms in Eq. 14 
we evaluate 


aho 


tc 

l + f c 


l f e,0 

& a o nD,o = — 

1 + F c o 




r c,N 

1 + f C ,N 

(18) 


Method 1 corrections to a baseline CFD solution 
are made using 


Qw{lnew(T)) = q w (CFD baseline) (19) 

(Method ljjnewiT)) - q w (Method 1,7 dd{T)) 


In the present application of Method 1, no attempt is 
made to account for changes to CFD extracted quanti- 
ties at the wall corresponding to the new surface catal- 
ysis model. 


Local, Discontinuous Change in Nose Region, Method 2 


More complete details of the following analysis 
may be found in a companion paper. 1 The post- 
processing algorithm used to correct for changes in sur- 
face heating associated with an abrupt, local change of 
the surface catalysis model in the nose region of a hy- 
personic vehicle is denoted Method 2. It is assumed 
that a change in TPS material occurs at a juncture 
(jump) location £j. Initially, it will also be assumed 
that the boundary-layer gas phase can be treated as 
frozen. Upstream and downstream of , both a e and 
r c vary rather slowly compared with a w itself in the 
relaxation zone, and hence are well approximated by 
local constants. The theory predicts a variation in oc w 
under these constraints as 

«* = cr< 1+f ^/ 2 + (£ > o) (20) 


Even though f c is discontinuous across the juncture 
the variation of a w across the juncture must be contin- 
uous (though its gradient will be discontinuous). Con- 
sequently, the constant C is proportional to the jump 
in r c and defined by 


C = 




(-U 

u + rr 


nW) 


,(l+f+)/2 


(21) 


Examination of Eqs. 20 and 21 reveals two features of 
the relaxation process. First, the further downstream 
the jump is, the longer it takes to adjust to the new 
catalytic condition. Second, the more catalytic the 
surface, the shorter the relaxation distance for any im- 
posed AT C . 

The diffusional contribution to the heat transfer 
function Q w will jump discontinuously across the junc- 
ture because of the discontinuous change in the gradient 
da w j dr] ^ r c^w • 


Qw =0.47PrV3[!_^ + ^i + Le -i/3^] 

QA7 d^l hD [ f ‘( Le " U - i] (22) 

The application of Eq. 22 in Method 2 to correct CFD 
baseline results obtained without the local change in 
surface catalysis follows the example set in Method 1. 
Note that in Eq. 22 the first term represents the con- 
tinuous part of the solution across the juncture and 
the second part represents the jump effect associated 
with the diffusional term. When the change in Q w is 
computed between the baseline case and the case with 
a jump in surface catalytic efficiency the contribution 
from the continuous part exactly cancels. 

Though not presented here, the assumption of 
frozen gas phase in the boundary layer can be elim- 
inated. In this case, a local iteration is required to 
account for the sudden change in boundary condition 
on the reaction integral. Inclusion of this effect was not 
significant in the test cases which follow. 
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Generalization for Complete Vehicle, Method 1G 

Equations 13-14 must retain the normalized ki- 
netic energy terms when treating flow away from the 
stagnation region. Equations 11-14 in Method 1 em- 
ploy the limiting form of the parameter \ the stag- 
nation point and utilize this result in the entire nose 
region. This approximation fails away from the nose 
region, and a better evaluation of parameter x ls re ~ 
quired. Numerical approximations to \ based on flat 
plate theory for downstream locations are also unsatis- 
factory. The available CFD solution may be exploited 
by solving for \ such that the integral boundary-layer 
solution for q w matches the known CFD solution for 
q w . This treatment avoids the numerical complexity of 
a detailed integration of Eq. 7 along streamlines. This 
same value of x is then used in the integral boundary- 
layer formulation for the new values of catalicity (and 
wall temperature, if required) to compute the new heat- 
ing distribution. The approximation for x ignores inte- 
grated effects of changes in surface temperature along 
the surface streamline, but works well in limited tests 
performed to date. 

In the most general case, viscous dissipation of ki- 
netic energy tends to produce a local maximum tem- 
perature within the boundary layer away from the stag- 
nation region. The higher temperature may lead to 
dissociation dominated chemistry within the boundary 
layer. Inclusion of this effect leads to the following cor- 
rection which is added to Eq. 14: 


A Q w , F = -QA7Pr l J 3 T D I D -^- 1 


where 7p is a reaction integral, discussed 
by Inger. 2 


Le w r c 
1 T f c 


(23) 

and tabulated 


Test Cases 

Two sets of tests for hypersonic flow over a 0.6 m 
radius sphere are presented. In the first set of tests, 
surface heating rates are computed for three differ- 
ent wall catalytic boundary conditions at trajectory 
points defined in Table 1. The wall temperature in each 
case is defined by a radiative equilibrium wall bound- 
ary condition computed with 7 = 7 (T) and emissivity 
6 = 0.9. Both the fully catalytic (7 = 1) and non- 
catalytic (7 = 0) solutions were obtained with the 
wall temperature distribution obtained in the finite- 
catalytic (7 = 7(T)) case. This specification was main- 
tained to eliminate variations in wall temperature as 
a cause of differences between the CFD and the inte- 
gral boundary-layer heating predictions. Actual varia- 
tions in radiative equilibrium surface temperatures that 


would be expected in the general case would be handled 
by a local iteration to obtain self consistent heating and 
temperature distributions. 


Table 1 - Trajectory and Heating 


Pt 

1/00, 

m/s 

poo , 

kg/m 3 

Too, 

K 

<h, 

5 

W/cnP 

?o, 

1 

5493. 

2.00 “ 4 

238. 

50.5 

41.1 

24.4 

2 

4440. 

3.53 

250. 

34.3 

30.2 

20.2 

3 

3551. 

6.64 ~ 4 

264. 

22.2 

18.6 

14.6 

4 

2974. 

1.01 - 3 

271. 

15.5 

14.6 

— - 


In the second set of tests, an abrupt change in sur- 
face catalysis (from finite-catalytic to fully-catalytic) is 
introduced on the sphere for trajectory point 1 of Table 
1 at s/R = 7t/6. Such a change would approximately 
correspond to use of a glassy coated surface at the 
stagnation point followed by polished (non-oxidized), 
metallic surface. The integral boundary-layer method- 
ology is used to predict both the overshoot in heating 
rate at the juncture and the relaxation distance to the 
new heating level corresponding to 7 = I. Wall tem- 
perature distributions are held fixed at the radiative 
equilibrium values corresponding to the finite-catalytic 
solution as described in the previous paragraph. 


Global Changes to Surface Catalysis 

In Fig. la, the non-catalytic solution is treated as 
the available CFD baseline solution. This baseline CFD 
data set has been post-processed using Method 1 to 
obtain corrected heating levels associated with changes 
in the surface catalysis model to finite catalytic and 
to fully catalytic. The actual CFD heating predictions 
with the new surface catalysis models are compared 
with the results of Method 1 in Fig. la. Comparisons 
are generally very good, with a maximum difference in 
predicted stagnation point heating level occurring for 
7 = 1 case equal to 10%. 

Similar comparisons are made in Fig. lb for point 
1 except that the finite-catalytic CFD solution is used 
as a baseline and post-processed corrections are derived 
for the fully-catalytic and non-catalytic surfaces. Com- 
parisons with CFD heating predictions are within ap- 
proximately 10% over the entire surface. The fully- 
catalytic CFD solution is used as a baseline in Fig. 
lc. The finite-catalytic prediction from Method 1 is 
within 10% of the CFD prediction using this baseline; 
however, the non-catalytic prediction is approximately 
30% low as compared to the CFD result. Alternative 
methods for extracting and averaging boundary-layer 
quantities to improve this comparison while maintaing 
good comparisons from other baselines across the entire 
trajectory continue to be explored. The Lewis number 
profile across the boundary layer is sensitive to the cat- 
alytic boundary condition and an additional correction 
within the theory may be advised. 


5 




(a) CFD baseline at 7 = 0 


(a) CFD baseline at 7 = 0 



(b) CFD baseline at 7 = 7 (T) (b) CFD baseline at 7 = 7 (T) 



(c) CFD baseline at 7 = 1 

Figure 1: Comparison of CFD heating levels (symbols) 
with analytic corrections (lines) derived from CFD 
baseline for point 1 as function of three surface catalysis 
models. 



(c) CFD baseline at 7 = 1 

Figure 2: Comparison of CFD heating levels (symbols) 
with analytic corrections (lines) derived from CFD 
baseline for point 2 as function of three surface catalysis 
models. 
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(a) CFD baseline at 7 = 0 


(a) CFD baseline at 7 = 7(X) 






(b) CFD baseline at 7 = 7(F) 


(c) CFD baseline at 7 = 1 

Figure 3: Comparison of CFD heating levels (symbols) 
with analytic corrections (lines) derived from CFD 
baseline for point 3 as function of three surface catalysis 
models. 


(b) CFD baseline at 7 = 1 

Figure 4: Comparison of CFD heating levels (symbols) 
with analytic corrections (lines) derived from CFD 
baseline for point 4 as function of three surface catalysis 
models. 

Comparisons of Method 1 with CFD at successive 
trajectory points are presented in Fig. 2 for point 2, 
Fig. 3 for point 3, and Fig. 4 for point 4. A non- 
catalytic baseline was not run for case 4 because effects 
of catalysis are so small at this point in the trajectory. 
In all of these tests, Method 1 corrected heating lev- 
els are judged to be acceptable for preliminary design 
purposes. 

Local, Discontinuous Changes to Surface Catalysis 

In Fig. 5, the finite-catalytic solution is treated as 
the available CFD baseline solution at trajectory point 
1. Method 2 is used to post-process the CFD solution 
to predict the heating levels associated with a local, 
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discontinuous change in surface catalysis to a fully cat- 
alytic condition at s/R = i r/6. This result is compared 
to a corresponding CFD computation with identical 
wall temperature distribution and grid density. The 
CFD results in Fig. 5a are computed with 30 cells in 
3 degree increments (A s/R = 7r/60). This grid density 
would generally be considered adequate to resolve the 
flowfield over a hemisphere in the streamwise direction. 
The integral bound ary- layer theory (Method 2) shows 
a spike in predicted heating level of thickness equal to 
one cell and magnitude more than double the upstream 
location. (The spike is somewhat rounded because the 
plotting routine averages cell centered quantities in or- 
der to plot at cell corners.) The CFD solution shows 
only a mild overshoot in heating on this grid. While 
a discontinuous rise in heating is expected across the 
interface, the relaxation process back down to the fully 
catalytic level is not properly resolved. The lack of 
resolution is not a limitation of the theory; rather, the 
theory is only applied at CFD grid points in the present 
algorithm. 

A factor of four increase in resolution is introduced 
in Fig. 5b. Even on this scale, the single point spike 
of Method 2 indicates that the relaxation process is 
occurring on a scale still smaller than the current grid 
resolution (A s/R = 7r/240) . The CFD solution shows 
a larger overshoot on the finer grid, but is still roughly 
half the value predicted by Method 2. 

An additional, local grid refinement providing a 
factor 10 more resolution in the vicinity of the inter- 
face is applied. The CFD solutions are compared on 
an expanded scale in Fig. 5c. A relaxation process is 
finally evident on this scale, though grid convergence 
has not yet been confirmed. The CFD heating over- 
shoot is approximately 25% larger on this refined grid, 
but still approximately 33% smaller than the Method 
2 prediction. Resolution of these differences is still be- 
ing worked; nevertheless, the application of Method 2 
clearly identified a grid resolution problem and serious 
underprediction of the heating overshoot of the original 
CFD solution, which rose from 50 to 80 W/cm 2 in the 
course of two grid refinements. 

The present test case provides an extreme exam- 
ple of the overshoot that can be obtained with abrupt 
changes in surface catalysis. Usually, the juncture 
would be placed much further from the stagnation 
point. Theory shows that the relaxation distance in- 
creases with running length from the stagnation point, 
easing somewhat the necessity for local grid refinement. 
Also, diffusion in a streamwise direction, across the dis- 
continuity, has not been included in the analysis but 
should serve to reduce the heating spike. Nevertheless, 
the integral-boundary-layer analysis provides a-priori 
guidance as to the spacial resolution required to com- 
putationally address these issues. 





(c) Local refinement with A s/R — n/2400. 

Figure 5: Comparison of CFD heating levels obtained 
for an abrupt change in surface catalysis from 7 = 7(T) 
to 7 = 1 with analytic corrections derived from CFD 
profiles obtained at 7 = 7(T) for point 1. 
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Access-to-Space Vehicle Application 



(a) CFD baseline at 7 = 0 



(b) CFD baseline at 7 = 7(T) 



Figure 6: Comparison of CFD heating levels (symbols) 
with Method 1G predictions (lines) for RLV reference 
sphere at three catalytic surface conditions. 


All previous test cases focused on the stagnation 
region of a sphere where the limiting form of param- 
eter x ls easily defined. In this section, application 
of Method 1G to simple and complex configurations is 
tested. The boundary-layer edge is defined as the point 
where total enthalpy equals 0.995 of total freestream 
enthalpy. Velocity at the boundary layer edge is taken 
as the component that is perpendicular to the surface 
normal. The Lewis number is calculated with a diffu- 
sivity of atomic oxygen diffusing in molecular nitrogen. 

In the first case, Method 1G is applied to analysis 
of heating to a sphere at Mach 13.5 and 58 km altitude. 
Heating rates are plotted as a function of streamwise 
computational coordinate i in Fig. 6. The CFD cases 
are fully converged on grids of 128 cells across the shock 
layer and 40 cells around the body. The CFD solution 
for each surface catalysis model was computed using 
the same wall temperature distribution, corresponding 
to radiative equilibrium with the fully catalytic sur- 
face. The integral-boundary-layer corrections to any 
CFD baseline using Method 1G are seen to be good pre- 
dictors of CFD heating levels as a function of surface 
catalysis away from the stagnation point (i > 20) as 
shown in Fig. 6. However, in the vicinity of the stagna- 
tion point the Method 1G predictions differ from CFD 
values by as much as 20% when using a non-catalytic 
CFD baseline. Predictions are much improved when 
using a catalytic CFD baseline. 

The second test case examines the front and mid- 
sections of a Reusable-Launch- Vehicle (RLV) candidate 
configuration for a flowfield simulation at Mach 25, 45 
deg. angle of attack, and 79.6 km. The surface grid 
representation for the front section is defined with 52 
streamise cells and 64 circumferential cells. The surface 
grid representation for the mid-section is defined with 
18 streamwise cells and 77 circumferential cells. The 
shock layer is resolved with 64 cells. Heating rates for 
both a fully catalytic wall and a finite-catalytic wall 
were computed at the respective radiative equilibrium 
wall temperatures. A secondary iteration was employed 
in Method 1G to accomodate the implicit, functional 
dependence of T c and 9 W on the wall temperature. 

Circumferential heating distributions as a function 
of computational coordinate j varying from leeside ( j — 
1) to windside ( j = 64) are presented in Fig. 7a for the 
front section. The circumferential cut is from the i = 
40 plane which lies far downstream from the nose and 
upstream of a wing. The windside centerline heating 
distribution as a function of computational coordinate 
i varying from the stagnation point (i = 16) to the exit 
plane of the front section (i = 52) is presented in Fig. 
7b. Method 1G corrections to CFD baseline heating 
are generally within 5% of computations. Discepancies 
in the stagnation region ( i < 20) are larger but still 
within 20% for the fully catalytic case. Here again, 
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(a) i = 40 circumferential distribution 



(b) j = 64 centerline distribution 

Figure 7: Comparison of CFD heating levels obtained 
on front section of RLV with analytic corrections de- 
rived from CFD profiles obtained at 7 = y(T) and 
7 = 1 at Mach 25 and 79.6 km. 

Method 1G corrections are in better agreement with 
CFD when computed using a fully catalytic-baseline. 

The mid-section of the RLV includes the windside 
surface, the wing, and a small part of the leeside sur- 
face. The CFD and integral-boundary-layer results for 
heating on this section are compared in Fig. 8. Com- 
parisons are generally within a few percent in a circum- 
ferential direction around the wing leading edge (Fig. 
8a), around the windside surface (Fig. 8b) and in a lon- 
gitudinal direction along the windside centerline (Fig. 
8c). Comparisons along the wing leading edge show 
differences that are generally less than 10%. 

These results are the first demonstration of the 
method for a complex flow. Even in the case where 
significant turning of streamlines occurs for flow ex- 





(c) j = 19, 77 centerline and wing leading edge 

Figure 8: Comparison of CFD heating levels obtained 
on mid-section of RLV with analytic corrections derived 
from CFD profiles obtained at 7 = 7(T) and 7 = 1 at 
Mach 25 and 79.6 km. 
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panding around the side of the vehicle (32 < j < 48) 
in Fig. 7a the integral-boundary-layer extrapolation 
from the baseline CFD computation is an excellent pre- 
dictor of the CFD result at the off-baseline catalysis 
model. Though not shown in the figures, the Method 
1G predictions for a non-catalytic surface using either 
the fully-catalytic baseline or the finite-catalytic base- 
line also agree within a few percent. The largest er- 
rors in Method 1G occur in the vicinity of a stagnation 
point when using a non-catalytic CFD baseline. Some 
preliminary perturbations of this model suggest that 
the manner in which nitrogen recombination is handled 
may be the error source. 

The results in Fig. 7 and 8 do not include the 
correction for dissociation dominated chemistry in the 
boundary layer as defined in Eq. 23. The correction 
was evaluated and found to represent less than 1% of 
the total heating along the windside centerline for this 
case. The term grew unreasonably large as the leeside 
was approached. It appears that the curve fits to tabu- 
lated data for the reaction integral 2 ip are used outside 
their range of validity. It is expected that this term 
would be of greater significance if the vehicle were at a 
smaller angle of attack and more dissociation occurred 
within the boundary layer. 

The differences shown for this test case are thought 
to be small enough to be acceptable for preliminary 
thermal protection system design purposes. Some ad- 
ditional tuning of the method is still required in the 
stagnation region. Also, the consistency of these pre- 
dictions over a broader range of entry conditions and 
geometric complexity remains to be established before 
Method 1G (or a related approach) can be used with 
confidence in a design mode. The calculations pre- 
sented here using both Method 1 and Method 1G offer 
strong evidence that this design goal is within reach. 


Concluding Remarks 

A new approach for combining the insight af- 
forded by integral boundary-layer analysis with com- 
prehensive (but time intensive) computational fluid 
dynamic (CFD) flowfield solutions of the thin-layer 
Navier-Stokes equations is presented. The approach 
extracts CFD derived quantities at the wall and at the 
boundary layer edge for inclusion in a post-processing 
boundary-layer analysis, valid for most of the acreage 
on a realistic hypersonic vehicle. It allows a designer at 
a workstation to address two questions, given a single , 
baseline CFD solution: (1) How much does the heating 
change for a thermal protection system with different 
catalytic properties than was used in the baseline CFD 
solution? (2) How does the heating change at the in- 
terface of two different TPS materials with an abrupt 


change in catalytic efficiency? 

Global changes in surface catalysis lead to global 
changes in heating rates. Two methods were developed 
in order to quantify these changes from a baseline, CFD 
solution. Method 1 is restricted to the stagnation re- 
gion of vehicles. Method 1G is a straightforward ex- 
tension of Method 1 which is applicable to the entire 
vehicle. Rather than employing limiting forms of the 
boundary-layer equations or integrating along stream- 
lines, it solves for an integral-boundary-layer streamline 
parameter to match the local CFD baseline heating 
rate. This parameter is then used to predict heating 
rates at off-baseline values of surface catalysis. The 
influence of boundary-layer chemistry can be included 
in these analyses, but was found to be small in the 
test cases treated herein. Off-baseline predictions of 
heating rates by both Method 1 and 1G were gener- 
ally in very good agreement with CFD solutions, with 
Method 1G appearing somewhat superior in the limited 
tests performed to date. The differences shown for this 
test case are thought to be small enough to be accept- 
able for preliminary thermal protection system design 
purposes. Some additional tuning of the method is still 
required in the stagnation region. Also, the consistency 
of these predictions over a broader range of entry con- 
ditions and geometric complexity remains to be estab- 
lished before Method 1G (or a related approach) can 
be used with confidence in a design mode. The calcula- 
tions presented here using both Method 1 and Method 
1G offer strong evidence that this design goal is within 
reach. 

Abrupt, discontinuous changes in surface cataly- 
sis lead to large, local changes in heating rates that 
must relax to the new level dictated by the catalytic 
efficiency of the new surface. Method 2, applicable 
only in the stagnation region, predicted large heat- 
ing spikes associated with an abrupt change in sur- 
face catalysis. The spikes were larger than initially 
predicted by CFD. However, the relaxation zone, as 
defined by the integral-boundary-layer method, was 
not adequately resolved and grid refinement tended to 
bring the CFD result in better agreement with that 
predicted by Method 2. Remaining differences between 
Method 2 and CFD are likely due to numerical diffusion 
across the discontinuity and require further study. An 
extension of Method 2 to Method 2G, applicable out- 
side the stagnation region, should be straightforward 
but has not yet been tested. 
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